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Abstract 



< 

The thermonuclear burn-up of highly compressed deuterated methane CD4 is considered in 



the spherical geometry. The minimal required values of the burn- up parameter x = po ■ rj are 

re- 
determined for various temperatures T and densities po- It is shown that thermonuclear burn- up in 

CD4 becomes possible in practice if its initial density po exceeds ~ 5-10 3 g-cm~ 3 . Burn-up in CD2T2 

a' 

methane requires significantly (~ 100 times) lower compressions. The developed approach can be 



used in order to compute the critical burn-up parameters in an arbitrary deuterium containing 
fuel. 
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In this communication the thermonuclear burn-up in highly compressed CD4 (and CD2T2) 
methane is considered. Our goal is to determine the minimal required values for the burn- 
up parameter x (where x = po ' r /)> density po an d temperature T at which thermonuclear 
burning starts locally and then propagates successfully to the rest of the highly compressed 
CD 4 methane. Let us assume that initially we have the infinite, homogeneous, immovable 
CD4 (and CD2T2) methane (gas), which is assumed to be cold at time t = 0. The burn- up 
problem is considered in the following form. Inside homogeneous immovable CD4 methane a 
quite large amount of energy Qo is instantaneously released in a very small (point) volume. 
The temperature in this volume increases very rapidly to extremely large values. As a 
result a high-temperature thermal wave forms, which begins to propagate from this very hot 
volume into the rest of the cold methane. After some time this wave transforms into the 
shock wave, which propagates and disappears, in the general case, in the cold methane. 

However, if the three following conditions are obeyed: (1) the initial density of the highly 
compressed methane po is quite large, (2) the temperature in the central, heated region 
exceeds the ignition temperature T > T c and (3) its volume is large enough V > V c , then 
the thermonuclear burn wave may form and propagate out, igniting the rest of the cold CD4. 
In general, the thermonuclear burn wave propagates initially as a thermal wave, and later 
as a detonation wave, i.e. the shock wave which supports itself 

In this sudy we want to determine the critical conditions, i.e. the minimal possible po,T c 
and V c values at which the combustion zone still expands to the infinity in the cold CD4 
or CD2T2 methane, respectively. In the general case, this problem is very complicated, but 
here we restricted ourselves to the spherical geometry only. This means that the initial hot 
spot has a spherical form, and later, the expanding thermal or detonation waves also remain 
spherically symmetric. In this case the arising gas flow is one-dimensional, and therefore, the 
appropriate equations can be also written in the one-dimensional form (see below). Finally, 
the initially complicated problem of V c determination is reduced to the one dimensional and 
relatively simple problem: to find a critical value for some linear r c parameter, which can be 
called by the critical radius of the initial hot spot. The consideration of the cylindrical and 
planar cases can be based on analogous assumptions, but our present study is concentrated 
mainly on the spherical symmetry, since this case is of specific interest for various applications 
in practice. 

As is well known (see e.g. high-temperature burn-up in thermonuclear fuel is 
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governed by the so-called burn-up equation. Let 77 be the space radius of the hot (i.e. 
combustion) zone at the time t. The burn-up parameter x = po ■ r/ (where po is the initial 
density) plays a very important role below. The initial values of rf and x are designated 
below as r = r c and Xq = x c . This corresponds to the physical sense of r c as the minimal 
or critical radius of the hot zone for which the thermonuclear burn wave still propagates 
to infinity (in the cold CD4 methane). The temperature T behind the shock or thermal 
wave (i.e. in the hot or combustion zone) is assumed to be significantly higher than in the 
initially cold CD4 methane, where To ~ 0. Furthermore, without loss of generality we may 
assume that T inside of the hot zone does not depend upon spatial coordinates. Likewise, 
the ion and electron temperatures are equal to each other and coincide with T. Such an 
approximation has a quite good accuracy for burn-up problems, unless the temperature T 
exceeds 25 keV In this approximation the burn- up equation can be written (for detail 
see e.g. [2]) in the form which contains only one variable, namely the burn-up parameter 
x = po ■ Tf (g • cm~ 2 ): 

dT = v T | q{x,T) 

dx X C • Vmax 

where q(x,T) is the so-called normalized energy release function, v — 3, 2 and 1 for the 
spherical, cylindrical and plane cases, respectively. C is the specific heat per unit mass, 
which is assumed to be a constant on the spatial coordinates and time. V mSLX = is the 
maximal velocity of the hot (i.e. combustion) zone expansion. Since such a zone expands 
either by the detonation or thermal waves, the velocity V mSiX is the larger of the two cor- 
responding velocities, i.e. V max = max{(^-)d, (~if)t}- In general, it can be shown that 
the thermal wave may propagate the thermonuclear burning only at hardly attainable tem- 
peratures T > 23 - 25 keV and only in the dense equimolar DT— mixture. In all other 
cases, including DT— mixtures with low (< 10%) tritium concentrations, highly compressed 
deuterium, DT— hydrides and deuterides of light elements, etc, the thermonuclear burning 
(at T < 25 keV) propagates by using only the relatively slow, high-temperature detonation 
expansion. This means that without loss of generality we can assume in the last equation 
Vmax = Vd, where Vd is the velocity of the detonation wave. 

To determine the velocity Vd we shall apply the strong explosion approximation, which 
has a good accuracy for the considered problems (see e.g. [2I). In this approximation 
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the expression for Vd takes the from 
V D = a( 7 ) • VCT~T 



3-7T- (7 + I) 2 • (7 



— • ■ T = 5.28443640 ■ 10~ 3 • VC ■ T (2) 



\ 8.(3-7-1) 

where 7 = | and C is in MJ • fir 1 • /ee^ 1 (IMJ = 1 ■ 10 6 J), T is in fceF and V R isjn 
cm ■ nsec~ x . The specific heat C (in MJ ■ g~ x ■ keV~ x ) per unit mass is of the form 
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].Q: 



c= 144.7164346 M J ■ ■ fceV -1 (3) 

where the A and Z values are the mean atomic mass number and the mean atomic number, 
respectively. They must be expressed in units in which the proton charge and proton mass 
equal 1. Also we assume that the difference between the proton and neutron masses is 
negligible and also the mass defect for all considered nuclei equals zero exactly. Both these 
approximations are traditional in thermonuclear problems, and they have been widely used 
in earlier works (see e.g. Q or fl). Therefore, for the considered CD 4 methane Z = 2.00 
and A = 4.00, while for the CD2T2 methane Z = 2.00 and A = 4.40. From here one easily 
finds: C = 108.537260 MJ ■ g~ l ■ keV' 1 , V D = 5.50539183 ■ 10~ 2 ■ VT cm ■ nscc~ x for CD 4 
methane and C = 98.6702963 MJ ■ cT 1 • keV' 1 , V D = 5.24918515 ■ 10~ 2 -VT cm- nscc~ l for 
CD2T2 methane (T is in keV). 

To discuss the energy release function q(x, T) in Eq.flTJ we introduce the useful pair nota- 
tion (a,f3) which designates the (d, <£)— pair for CD4 methane and (d, t)— pair for CD2T2 
methane. By using this notation the corresponding energy release function q(x, T) = 
q a< p(x,T) Eq.(|T]) can be written in the following universal form (for more detail see [2J): 

q a> p(x, T) = [i p ■ s a> p(x, T) ■ (1 + B a> p ■ — — ) - H (4) 

1 + o a ,px 1 + c Qj/3v /p • xT 4 

where [i at p is the ratio of the (a, /3) fragment mass to the total atomic (quasi-molecular) 



mass and the factor K a> p takes the well known form jg]: 

1 (Zinj-ZiKZimZ?) 1 (Eini-Z^XiTiiZf) 
Ka '" Kp(J2 a n a Z a )& a n a Zl) Nl, 2 2 

where the summation is taken over all ions present. The factor j^— = 7J2 = \ corresponds 

a, /3 

to the fact that CD4 contains 4 deuterium atoms or two d — d fragments (i.e. Ndd — 2). 
Analogously, there are two d — t fragments in each molecule of CD 2 T 2 methane (i.e. = 
2). For instance, from the last equations for the CD 2 T 2 methane one easily finds that 



i/ — n — 2+2+3+3 _ 10 QT .J ., _ _ 1 6+1+1+1+1 G^^ + l^ + ^ + l 2 _ 9r . p . „ 
/^*,/3 — /M* — 12+2+2+3+3 ~~ 22 aI1Cl 'V- 8 ~~ ~~ 4 ' 2 ' 2 ~~ ^ mall Y: 

K d,t ■ Hdt = 5.165289256. The use of these two factors fi a ^ and K a ^ allows us to apply the 
same short (i.e. ionic) energy release functions s at p(x,T) and bremsstrahlung constants the 
C a fi factor as for the appropriate pure deuterium or equimolar DT-mixture. At least in the 
first order such an approximation seems to be quite good js]. In general, the parameters 
B a ,/3,b a ,i3 and c a ^ depend upon the ionic contents of the considered mixture. But for our 
present approximate purposes we shall assume that in the CD4 and CD2T2 methane they 
are exactly the same as for pure deuterium and equimolar DT-mixture, respectively. Their 
numerical values can be found in l|, 

The short energy release functions s ai p(x,T) in Eq. (Tj0) correspond to the energy gain 
produced only by positive high-energy ions arising in the thermonuclear reactions. They 
can be expressed also in the following general form: 

, MJ-cm 3 

Sa,p = Y ■ ((Tap ■ V a p) — (6) 

g z ■ nsec 

where, e.g. cr^t is the (d, t)— reaction cross-section, Vdt is the relative velocity of the deu- 
terium and tritium nuclei and (WJ t ■ Vdt) means the averaged (dimensionless) value over the 
whole range of relative velocities. The numerical factors Y are [21]: 8.17611 ■ 10 18 for the 
(d, t)— reaction and 3.03891 • 10 19 for the (d, d)— reaction when all reactions for 3 He nuclei 
are ignored. However, if the (d, d)— reaction is considered in the presence of high intensity 
neutron fluxes, then the factor Y equals 4.63191 • 10 19 . The last case represents the situation 
when all 3 He nuclei react with neutrons, and the following (d, t)— reactions are included in 
the consideration. Now, the remaining problem is to find an analytical or numerical ex- 
pression for the (add 4 v dd) an d (cr<& ' v dt) values which depend upon the temperature T only. 
Presently, they have been chosen from [2], but it should be mentioned that the analogous 
(&dt ■ Vdt) function from [1] and (add ■ i>dd) function from 7] produce quite close results. 

Now, the solution of the burn- up equation Eq.(JTJ not contain anyprinciple difficulties. 
Moreover, this solution can be written in the following integral form [2(: 



T «= T °'<-> V D ( v ,T(y)) < 7 > 

where xq = po ■ r o and T(xq) = To are the initial value of the burn-up parameter and temper- 
ature, respectively. The general T(x) dependence has the following qualitative behaviour, 
when x grows: it starts at x = x (T(x ) = T ), then decreases, reaches a minimum and later 



increases to infinity. The minimal temperature T m i n can not be negative and the appropriate 
condition T min = is used to determine the critical value xo = x c {= p ■ r c ) of the burn-up 
parameter. 

The results of our numerical calculations are presented in Table I (for CD4 methane) and 
Table II (for CD 2 T 2 methane). These Tables contain the x c values for various temperatures 
T (4 keV < T < 20 keV) and different initial densities po- In Table I the two subcases are 
considered for each temperature and density. The left subcolumn (a) in Table I represents 
the results when all reactions with 3 He nuclei are ignored. The right subcolumn (b) in this 
Table corresponds to the case, when all arising 3 He nuclei react with neutrons, i.e. the 
thermonuclear burn-up in the presence of high intense neutron fluxes. In the last case we 
assume that the (n, 12 C; 13 C, 7) reaction has relatively small cross-sections. As it follows from 
Table I the propagation of the thermonuclear detonation wave in CD4 methane simplifies 
significantly in the presence of high intense neutron fluxes. This effect is of paramount 
importance for relatively small (or realistic) temperatures T m 4 - 6 keV. In these cases the 
critical burn-up parameter drops by approximately 3 times. The thermonuclear detonation 
in CD 2 T 2 methane (see Table II) is not very sensitive to the presence (or absence) of any 
neutron fluxes. 

In discussing the results from Tables I and II, we note that the thermonuclear burn- 
up in CD 4 methane requires significantly greater compression than in CD 2 T 2 methane. For 
realistic temperatures T ~ 4 — 6 keV the deuterated CD4 methane must be compressed ~ 75 
- 100 times more strongly than CD 2 T 2 . In other words, the thermonuclear detonation in CD4 
methane may propagate successfully only if its initial densities po ~ 2.5 • 10 3 — 5 ■ 10 3 g ■ cm" 3 
or higher, which are quite close to the corresponding Fermi limits (see below). Analogous 
'critical' densities for CD 2 T 2 methane are ~ 70 - 90 g ■ cm" 3 and quite comparable with 
the densities required in order to produce thermonuclear burn-up in the dense deuterium or 
DT— mixtures with very low (< 1%) tritium concentrations (see 2|). It should be mentioned 
also that the minimal energy £ c , which is needed to produce thermonuclear ignition in the 
spherically symmetric case, can be easily evaluated in terms of the known x c parameter or 
critical radius r c : £ c = 4^ ■ x 3 ■ po 2 ■ C ■ T = 4^ ■ r 3 • p ■ C ■ T, where po is the initial density 
and T is the burn-up temperature. 

Note also, that the extremely violated energy balance in the deuterated CD4 methane 
produces another complication which can not be found in CD 2 T 2 . The source of such a 
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complication is obvious. Indeed, the velocity of the detonation wave Vd increases with the 



significantly longer time than in CD 2 T 2 . Therefore, we would expect that at some critical 
temperature T cr the shock wave breaks away from the burn wave and propagates into the 
cold fuel without producing any ignition. Since the energy losses related with the expansion 
of the hot zone are not compensated, the temperature behind the shock wave decreases 
rapidly T ~ rj 3 . Later, when the spatial radius of the hot zone increases to the r cr value 
and its temperature drops to the T = T cr value, then the thermonuclear burn-up begins, i.e. 
the initial, very fast shock wave slows down and transforms into the detonation wave with a 
lower temperature. In CD 4 methane such a critical temperature T cr 11.5 — 13.75keV and 
its value depends on the initial density p . The existence of the critical temperature means 
also that in CD4 methane the critical burn-up parameter x c (T,po) has the finite, non-zero 
limit at T — > +00: lim*r_;. +00 x(T, p ) = x(T cr ,p ). In other words, for CD 4 methane the 
critical burn- up radius r cr {r cr = x ( Tcr ^ Po > ) can be determined uniformly for each density, and 
moreover, its value: (1) does not vanish when temperature grows, and (2) does not depend 
on the temperature T, if T > T cr . In fact, this means that the non-physical part of the 
burn-up curve (where > 0) is ignored and replaced by the straight line ^ = 0. For the 
CD 2 T 2 methane ^ < everywhere (at T < 20keV) and the principal difference between 
the thermonuclear burn-up in CD 4 and CD 2 T 2 can be easily understood from comparison of 
Tables I and II. The x and T cr values are also presented (for each density) in Table I. Since 
T cr < 13.75keV in all cases, the higher temperatures (T > 14 keV) are not considered in 



Now, it is necessary to show that the obtained solution corresponds to the actual burn-up 
process. In general, successful burn-up means that the two following conditions are obeyed: 
firstly, the energy release in the hot zone behind the expanding thermonuclear burn wave 
must exceed all energy losses related with the hot zone expansion. This gives the burn-up 
equation Eq.flTJ) used above. Secondly, the thermonuclear burn wave must move faster than 
the following disassembly (or rarefaction) wave. The later is essentially the high-temperature 
sound wave. Its velocity C s takes the following form: 



temperature as Vd — 0, • yT, while the thermonuclear ignition in CD4 takes obviously a 



Table I. 





where T is in keV, 7 



I and C is determined above. Pi 



Hi ■ T and P e 



n e ■ T are 
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the ion and electron pressures, respectively. As follows from Eq.Q in the strong explosion 
approximation we always have Vd(T) > C S (T) for arbitrary T. In other words, in this 
approximation the disassembly wave always moves behind the thermonuclear burn wave, 
i.e. in the combustion zone. The motion of the disassembly wave in the combustion area 
does not mean the termination of thermonuclear burning. However, if the disassembly 
wave passes the thermonuclear burn wave, then the fuel will be disassembled before the 
thermonuclear burning may start effectively. Such a situation can be found for both relatively 
small (~ IkeV) and very high (> lOOkeV) temperatures. 

In particular, at small temperatures and high compressions the electrons are very close 
to the Fermi limit of electron degeneracy. As a result, in this case, the velocity of the 
disassembly wave is significantly larger than it follows from Eq.flS]). Indeed, in this limit 
the electron pressure is given by the same formula P e = n e ■ T e , where the effective electron 
temperature T e takes the form: 

7T 2 ,T P ^ 7vr 4 ,T f 



T P = T, 



ef ■ 



1 + — •( — ) 2 + — •( — ) 4 + ... 

15 l T e / 150 l T e / 



(9) 



2 , x 2 



where T e f is the equivalent Fermi temperature, given by (see e.g. |8j): 

T ef (keV) = \ ■ fel • f- ■ = 9.3019938991 • lO" 3 • (|)S • ( Po )f (10) 

5 8tv 2m e A A 

where Na, h and m e are the Avogadro number, Planck constant and electron mass, respec- 
tively [9|. For instance, for 12 CD^ methane with p = 1 • 10 6 g ■ cm' 3 one easily finds 
T e f = 58.50 keV. Therefore, for T e = 1 keV the effective electron temperature T e m T e f 
and actual electron pressure P e exceeds its classical value in 58.60 times. In this case, the 
disassembly wave moves in w 5 times faster than it follows from Eq.flH]). Briefly, we can say 
that the appropriate correction on the electronic degeneracy is obviously needed when the 
actual temperature T(= T e ) is comparable with the equivalent Fermi temperature T e f ^]. 
In our present analysis the largest T e f value can be found in CD 4 methane, where po < 1 • 10 4 
g ■ cm -3 , i.e. T e f < 0.27199195 keV. This indicates that the maximal T e f temperature is 
significantly smaller than temperatures used in the present study (T > 4 keV), i.e. the 
appropriate Fermi correction to the electron pressure P e is very small and can be ignored. 

In conclusion, it should be mentioned that the presented approach predicts that the ther- 
monuclear burn-up is quite likely also in the higher deuterated alkanes (C n D 2ri+ 2), alkenes 
(C n D 2n ), alkynes (C n D 2n _2), etc. Indeed, the both factors K^a an d Hd,d have the finite limits 
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when n —¥ oo (these limits equal \ and ||, respectively). The determination of the critical 
parameters for each of these polymers is straightforward. Moreover, the critical burn-up 
parameters for an arbitrary deuterium containing fuel can be computed by using the formu- 
las presented above. In particular, for all deuterides of light elements with Z < 9 (D 19 F is 
included) such an analysis can be found in our paper [10]. 
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TABLE I: The critical values of the burn-up parameter x c (in g ■ cm~ 2 ) for various densities po 
(in g • cm -3 ), temperatures T (in keV) for CD4 methane. The superscripts mean that: (a) all 
reactions with the 3 He nuclei are ignored, and (b) the thermonuclear burn-up proceeds in very 
intense neutron fluxes. 





(a) 

Po 


(b) (a) (b) 

Po Po Po 


(a) 

Po 


o [h) 
Po 


T 


1 • 10 3 


1 • 10 3 5 • 10 3 5 • 10 3 


1 • 10 4 


1 • 10 4 


4 


13125.467 4330.014 2834.009 986.851 1542.655 564.471 


5 


9340.429 


3074.378 2019.777 703.628 1100.985 403.783 


6 


7783.949 


2552.553 1679.783 583.248 


913.755 


334.134 


7 


7010.860 


2288.540 1508.037 520.944 


817.637 


297.394 


8 


6585.202 


2138.550 1411.503 484.666 


762.567 


275.570 


9 


6341.346 


2047.877 1354.585 462.082 


729.256 


261.671 


10 


6206.550 


1992.443 1321.531 447.699 


709.104 


252.551 


11 


6145.359 


1960.466 1304.593 438.793 


697.844 


246.629 


12 


6135.744 


1945.657 1299.618 433.884 


693.165 


243.037 


13 


6135.744 


1943.435 1299.618 432.106 


692.812 


241.259 


X 


6135.744 


1943.435 1299.618 432.073 


692.812 


240.914 




11.62 


12.60 11.99 13.16 


12.38 


13.71 



TABLE II: The critical values of the burn-up parameter x c (in g ■ cm 2 ) for various densities po 
(in g ■ cm" 3 ), temperatures T (in keV) for CD2T2 methane. 



Po 


Po 


Po 


Po 


Po 


Po 


T 50 


100 


200 


T 50 


100 


200 



4 96.6663 60.2729 39.7280 10 8.4305 7.2564 6.1622 

5 50.6829 34.2076 24.0097 11 6.6828 5.9122 5.1495 

6 31.1823 22.4731 16.5634 13 4.5962 4.2287 3.8281 

7 20.9295 15.9492 12.2564 15 3.5037 3.3001 3.0624 

8 14.8525 11.8714 9.4686 18 2.6749 2.5653 2.4291 

9 10.9912 9.1474 7.5431 20 2.3703 2.2870 2.1811 
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